/)D~ /&<£& 


ANALYSIS  OF  MIXED  MODE  CRACK  GROWTH  IN  A  SHIP  HULL  MATERIAL 
UNDER  DYNAMIC  LOADING 


C.  R.  BARNES 
J.  AHMAD 
M.  F.  KANNINEN 

Stress  Analysis  and  Fracture  Section 
Battelle,  Columbus  Laboratories 

\  ABSTRACT 

^Research  on  fracture  In  HY-100  steel  under  mixed  mode  (Modes  I  and  II), 
dynamic,  elastic-plastic  loading  conditions  is  described.  A  combined  experi¬ 
mental  and  finite  element  analysis  effort  was  pursued.  The  analysis  employed 
the  J  criterion  for  fracture;  a  criterion  which  effectively  generalizes  the 
conventional  J-integral  to  mixed  mode  and  dynamic  loading  conditions.  The 
experiments  were  performed  with  inclined-notch  three-point-bend  specimens  of 
HY-100  steel.  By  coupling  the  experiments  and  finite  element  analyses,  a 

fracture  criterion  based  on  a  critical  lvalue  was  devised  and  verified. N 

i' 

INTRODUCTION 

The  fracture  mechanics  techniques  needed  for  materials  that  fracture  in 
a  highly  ductile  manner  must  give  explicit  attention  to  the  extensive  plastic 
deformation  surrounding  the  crack  tip.  Extensive  crack  tip  plasticity  has 
two  effects.  First,  because  of  crack  tip  blunting,  crack  growth  initiation 
requires  greater  loads  than  are  predicted  by  linear  elastic  fracture  mechan¬ 
ics  analyses.  Second,  significant  amounts  of  stable  crack  growth  under 
rising  load  can  precede  fracture  instability.  To  cope  with  these  essentially 
inelastic  processes,  elastic-plastic  fracture  mechanics  techniques  are  now 
being  widely  pursued  (1).  This  research  builds  upon  these  developments  to 
provide  the  basis  for  fracture  mechanics  treatments  under  the  more  general 
conditions  of  concern  to  ship  structures. 

Current  work  in  elastic-plastic  fracture  mechanics  is  largely  confined 
to  quasi-static  loading  conditions  and  to  crack  growth  in  the  opening  mode 
(Mode  I).  Yet,  in  an  impact  loaded  ship  structure,  crack  initiation  could 
occur  under  dynamic  conditions  with  a  combination  of  both  the  opening  mode 
and  the  sliding  mode  (Mode  II).  Thus,  it  is  necessary  to  obtain  a  quantita¬ 
tive  understanding  of  crack  growth  initiation,  propagation  and  arrest  under 
Initial  combined  mode  dynamic  loading  conditions.  The  research  described 
here  la  aimed  at  developing  this  understanding  through  a  program  of  inte¬ 
grated  experimental  and  finite  element  analysis  work.  The  objective  is  to 
provide  the  basis  for  the  more  general  inelastic-dynamic  fracture  mechanics 
approach  that  is  needed  for  Navy  applications. 

THE  GENERAL  APPROACH 

In  a  critical  survey  of  progress  on  elastic-plastic  fracture  mechanics, 
Kannlnen,  et  al  (1)  found  that  most  efforts  are  focused  on  the  J-reslstance 
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curve  approach.  Despite  the  fact  that  such  approaches  are  inherently  limited 
to  small  amounts  of  crack  growth  prior  to  fracture  instability,  they  do  pro¬ 
vide  a  significant  improvement  over  the  conventional  linear  elastic  fracture 
mechanics  techniques.  However,  the  current  approaches  are  restricted  to  Mode 
I  conditions  under  quasi-static  loading.  While  a  generalization  of  the  J- 
integral  valid  for  mixed  mode  and  dynamic  loading  conditions  exists,  it  has 
not  previously  been  critically  examined.  Nevertheless,  in  view  of  the  suc¬ 
cess  that  has  been  achieved  with  J  and  the  complete  lack  of  a  viable  alterna¬ 
tive,  this  parameter,  known  as  the  J-integral,  was  selected  for  use  in  this 

study.  The  mathematical  basis  for  J  and  the  manner  in  which  it  can  be 
applied  to  elastoplastic  dynamic  fracture  problems  can  be  found  in  references 
(2J ,  (3] ,  and  [4] . 

The  experimental  portion  of  the  research  described  in  this  paper  em¬ 
ployed  HY-100  steel  three-point -bend  specimens  with  slant  edge  cracks.  The 
specimen  configuration  is  shown  in  Figure  1.  This  configuration  has  a  number 
of  advantages  for  this  research.  First,  crack  Initiation  and  propagation  can 
be  obtained  for  both  quasi-static  and  dynamic  loading  in  the  same  specimen 
geometry.  Second,  by  varying  the  crack  angle  and  the  crack  length,  a  range 

of  Mode  I  and  Mode  II  combinations  can  be  achieved.  Third,  this  specimen  can 

be  economically  analyzed  with  a  finite  element  method.  Fourth,  because  of 
the  constraint  that  arises  in  bending,  the  J-reslstance  curves  obtained  may 
be  lower  bound  values  that  will  provide  conservative  predictions  when 
utilized  for  structural  integrity  assessments  [l]. 
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FIGURE  l.  THREK-POINT-BEND  SPECIMEN  WITH  A  SLANT  EDGE  CRACK 


A  dynamic  elastic-plastic  finite  element  analysis  was  used  in  the  anal¬ 
ysis  portion  of  the  research.  The  objective  was  to  learn  if  crack  growth 
initiation  in  HY-100  steel  is  governed  by  critical  values  of  J.  The  J  values 
were  determined  by  finite  element  analyses  of  the  experiments.  In  essence, 
by  calculating  the  vectorial  components  of  S  (i.e.,  and  Jj)  *•  *  func¬ 
tion  of  load,  and  using  the  load  level  at  initiation,  critical  values  for 
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experiments  performed  with  each  different  crack  angle  were  determined.  This 
provided  a  tentative  elastic-plastic  mixed  mode  fracture  criterion  in  terms 
of  a  single  material  parameter. 

Analyses  were  also  performed  for  impact  loading.  Here,  the  fracture 
criterion  deduced  under  quasi-static  loading  was  used  as  input  to  a  finite 
element  simulation  of  the  impact  loaded  experiment.  In  this  way  a  prediction 
of  the  Initiation  time  was  made  for  comparison  with  an  experiment.  The  quite 
reasonable  accuracy  of  this  prediction  then  gave  an  indication  of  the  useful¬ 
ness  of  the  approach  for  the  general  conditions  of  Interest  in  this  research. 

THE  EXPERIMENTAL  RESEARCH 


Quasi-Static  Loading  Experiments 

The  specimen  type  used  in  this  research  was  shown  in  Figure  1.  Speci¬ 
men  blanks  were  cut  from  HY-100  steel  plate.  The  chemical  analysis  and 
mechanical  properties  are  shown  in  Tables  1  and  2.  The  specimen  blanks  were 
oriented  such  that  the  crack  propagation  direction  for  the  specimens  with 
6  -  90°  (see  Figure  1)  was  parallel  to  the  plate  rolling  direction.  Crack 
propagation  in  specimens  with  6  -  30°  and  0  *  43°  was  at  60*  and  45°,  respec¬ 
tively,  to  the  rolling  direction. 

The  general  specimen  configuration  conformed  with  the  ANSI/ASTM  E604-77 
Standard  for  Dynamic  Tear  Specimens  except  for  the  thickness  and  notch  treat¬ 
ment.  While  the  ANSI/ASTM  standard  [5]  requires  a  thickness  of  0.623  inch, 
l.  15-inch  thick  specimens  (full  plate  thickness)  were  used  in  this  investiga¬ 
tion  to  promote  plane  strain  behavior  at  the  notch  tip.  The  standard  also 
requires  that  the  notch  root  be  sharpened  by  pressing  a  preformed  tool  into 
the  notch  root.  However,  past  experience  at  Battel le  has  demonstrated  that 
many  otherwise  ductile  materials  fail  either  partially  or  totally  by  cleavage 
when  specimens  containing  the  ANSI/ASTM  specified  notch  are  Impact  tested. 

TABLE  1.  HY-100  CHEMICAL  COMPOSITION 

(Average  of  Three  Samples) 
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TABLE  2.  HY-100  MECHANICAL  PROPERTIES 


Specimen 

Orientation 

2X  Offset 
Yield 
Strength 
(ksi) 

Ultimate 

Tensile 

Strength 

(ksi) 

Elongation 

(X) 

Reduction 
of  Area 
(X) 

Modulus 
(10^  psi) 

Longitudinal 

107.6 

121.3 

20 

72.4 

31.8 

(T-L) 

(742Mpa) 

( 836Mpa) 

(219Gpa) 

Trans ve rse 

108.1 

121.8 

20 

70.3 

31.1 

(L-T) 

(745Mpa) 

(840Mpa) 

(214Cpa) 

(One  such  material  Is  HY-100  steel).  In  addition,  the  specimen  configuration 
precluded  the  use  of  fatigue  to  sharpen  the  notch  root  (as  used  in  the  ASTM 
standard  for  K^c  specimens).  Consequently,  the  notch  was  sharpened  using  a 
jeweler's  saw  that  produced  a  root  radius  of  about  3.25  mils.  It  was  antici¬ 
pated  that  this  relatively  blunt  notch  would  store  enough  energy  to  promote 
rapid  unstable  crack  propagation  in  the  impact  tests  that  are  described 
later. 


Table  3  lists  the  specimens  tested  under  quasi-static  loading  condi¬ 
tions.  Tests  were  performed  using  a  Hat 'iwin-Taie  Emery  Test  Machine  having  a 
load  capacity  of  60,000  lbs.  The  loading  rate  was  adjusted  such  that  the 
load-point  displacement  (LPD)  was  0.01  inch  min~l.  The  loading  was  inter¬ 
rupted  at  preselected  points  on  the  load-displacement  curve  to  fill  the  notch 
with  a  alllcone  rubber  compound.  This  procedure  was  used  to  determine  the 
load  at  crack  initiation.  It  was  determined  that  the  cracks  initiated 
internally  and  tunnelled  for  a  distance  roughly  equal  to  the  shear  lip  width 
(0.2  to  0.3  inch)  before  reaching  the  surface. 


TABLE  3.  QUASI-STATIC  LOADING  EXPERIMENTS 


Specimen  Crack  Maximum 

Number  Position  6  Load  (lba) 


102S 

Mid-Point 

90* 

39,000 

(173.3  kN) 

101S 

45* 

42,500 

(189  kN) 

103S 

m 

30* 

43,000 

(173.9  kN) 
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Specimen  10 1 S  (B  *  45°)  was  farther  Instrumented  to  determine  strain 
versus  load  by  means  of  a  strain  gage  rosette,  A  crude  verification  of  the 
strains  was  provided  by  a  coarse  grid  applied  to  the  specimen’s  opposite 
face.  A  photographic  record  of  this  grid  was  made  at  each  load-strain  incre¬ 
ment.  An  identical  grid  was  also  applied  to  Specimen  103S  (B  *  30°)  and 
photographically  recorded.  Figure  2  shows  a  montage  of  Specimen  103S  from 
zero  load  to  fracture.  Figure  3  shows  the  fracture  surfaces  of  all  three 
specimens  tested  under  quasi-static  loading. 
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FIGURE  2.  SPECIMEN  103S  FROM  ZERO  LOAD  TO  FRACTURE 

The  experimental  data  from  tlic  three  quas l-stat leal ly  loaded  specimens 
Indicate  the  following.  First,  a  plot  of  load  versus  load-point  displacement 
(LPD)  shows  good  agreement  from  specimen  to  specimen  (Figure  4).  The  plot 
also  indicates  that  lower  values  of  B  have  higher  values  of  the  LPD  at  crack 
initiation.  Second,  the  crack-tip-opening  displacement  (CTOD)  at  crack  Ini- 
r{mr\nn  was  comparable  in  all  specimens  even  though  considerable  data  scatter 
was  evident  at  lower  loads  (Figure  5),  The  source  of  this  scatter  was  possi¬ 
bly  caused  by  deformation  zones  ooserved  about  4  mm  from  the  crack  tip. 
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FIGURE  3.  FRACTURE  SURFACES  OF  THE  QUASI-STATICALLY  LOADED  SPECIMENS 
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FIGURE  4,  LOAD  AS  A  FUNCTION  OF  LOAD-POINT-DiSPLACKHENT 
FOR  THE  QUASI-STATIC  LOADING  EXPERIMENTS 
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FIGURE  5*  LOAD  AS  A  FUNCTION  OF  CRACK  TIP  OPENING  DISPLACE  KENT 
FOR  THE  QUASI-STATIC  LOADING  EXPERIMENTS 


Dynamic  Loading,  Experiment! 

A  Battel le  designed  and  built  machine,  similar  in  concept  to  a  Charpy 
machine,  was  used  to  conduct  the  impact  tests.  It  is  a  multipurpose 
pendulum-type  aachine  with  a  total  impact  energy  of  approximately  16,000  ft- 
Ibs.  A  variety  of  specimen  geometries  can  be  accommodated  by  changing  the 
bolt-on  supports.  Figure  6  shows  a  photograph  of  the  specimen  loading 
arrangement  used  in  the  tests. 

To  assess  the  practicality  of  determining  crack  initiation  and  crack 
propagation  speeds  in  HY-100  steels,  a  crack  monitoring  gage  was  deposited 
upon  a  trial  specimen  surface  according  to  the  procedure  given  by  Kannlnen, 
et  ei  (6).  This  technique  Involves  the  vapor  deposition  of  an  electrically 
conducting  grid  upon  a  thin  epoxy  layer  bonded  to  the  test  specimen  surface. 

A  strain  gage  was  also  applied  to  the  specimen  at  the  locations  shown  In 
Figure  7.  In  addltton,  the  Impact  machine  was  electrically  wired  to  Indicate 
the  time  of  the  tup-specimen  contact.  The  tup  contact  switch  and  output  from 
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the  crack  monitoring  gage  were  recorded  simultaneously  by  a  Blomatlon  8100 
transient  recorder  and  an  FM  tape  recorder.  The  strain  gage  outputs  were 
recorded  on  two  adjacent  tape  channels. 

Although  the  data  from  this  feasibility  test  Indicated  that  the  dynamic 
Initiation  could  be  deduced  (Figure  8),  the  crack  speed  record  was  less 
definitive.  This  was  probably  caused  by  the  large  plastic  tone  generated  at 
the  crack  tip  which  resulted  in  epoxy-specimen  separation.  This  condition 
rendered  the  speed  measurement  suspect  since  past  experience  has  shown  that 
It  Is  essential  that  no  epoxy-specimen  debonding  occurs  for  reliable  crack 
speed  dete ruination.  This,  coupled  with  the  suspicion  that  crack  tunneling 
posed  similar  problems  In  determining  the  Instant  of  crack  initiation, 
prompted  specimen  modification  for  the  Impact  tests.  Accordingly,  three 
specimens  were  prepared  (one  each  with  8  ■  90* ,  45*  and  30*)  In  which  the 
notch  configuration  was  altered  as  shown  in  Figure  9. 


Tim#,  «wer©t#C0ft0l 


FIGURE  8.  RESULTS  OF  IMPACT  LOADED  FEASIBILITY  EXPERIMENT 
WITH  8  •  45*  (Tims  0  •  tup  contact) 
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Strain  gaga 
at  notch  center 


PICURE  9.  IMPACT  SPECIMENS  WITH  CRACK  STARTER 
NOTCH  STRAIN  CAGE  APPLICATION 

The  notch  geometry  shotm  In  Figure  9  vu  adopted  to  that  a  strain  gage 
could  be  applied  In  the  notch  center.  This  Indicates  crack  Initiation  sore 
precisely  than  would  a  gage  placed  upon  the  specimen  surface.  Another  gage 
vas  positioned  at  the  side  of  the  sharpened  notch  as  In  previous  specimens. 
This  gage  was  used  to  detect  the  time  at  which  the  Internally  Initiated  crack 
surfaced.  As  will  be  shown  later,  the  crack  Initiation  data  agreed  very  well 
with  the  data  fro*  the  feasibility  test,  which  Indicated  that  the  crack 
surfaced  at  about  200  uaec.  The  surface  crack  Initiation  times  for  6  *  90*, 
45*  and  30*  were  230,  222,  and  250  usee,  respectively.  Note,  however,  that 
the  feasibility  test  was  conducted  on  s  thinner  specimen.  Table  4  lists  the 
three  test  specimens  and  the  energy  required  to  fracture  the*. 

TAftLE  4.  IMPACT  LOADED  EXPERIMENTS 


Spe  clean 
Number 

Crack 

Position 

6 

Impact  Energy  (ft-lbs) 

1120 

Mid-Point 

90* 

1479  (2006  J) 

1 13D 

Mid-Pot nc 

45* 

1982  (2688  J) 

M4D 

Mid-Point 

30* 

1845  (2502  J) 

THE  ANALYSIS  RESEARCH 


The  analysis  effort  In  this  research  involved  performing  elastopUst lc , 
quasl-statlc  and  dynaalc  finite  eleaent  analyses  of  the  experl*ents.  The 
quasl-statlc  analyses  were  performed  to  coepute  the  1  Integral  values  corre¬ 
sponding  to  the  experimentally  determined  loads  at  crack  growth  Initiation, 
Using  the  J  values  at  Initiation  under  quasl-statlc  loading  conditions. 
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a  crack  Initiation  criterion  was  devised  and  used  in  the  dynaalc  analyses  for 
the  prediction  of  crack  Initiation  In  lapact  loading. 

All  finite  element  analyses  were  perforaed  using  Battelle's  elasto- 
plastlc  dynaalc  finite  eleaent  fracture  analysis  computer  code.  The  finite 
eleaent  aeshes  used  In  the  analyses  are  shown  In  Figure  10.  The  nuabered 
dotted  lines  In  Figure  10  indicate  the  three  contours  used  for  computing  the 
J-lntegrals.  The  values  obtained  using  the  three  contours  were  found  to  be 
fairly  path  Independent  [A].  The  results  of  the  quasl-statlc  analyses  sum¬ 
marised  in  Table  5  represent  the  average  of  the  three  contours  at  Initiation 
load.  K 


FIGURE  10.  FINITE  ELEMENT  MESH  FOR  30\  4S#,  and  90 #  SPECIMENS 
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TABLE  5.  ?!  and  J2  AT  INITIATION  FOR  QUASI-STATIC  LOADING 


■M  •  ««1 

a 

?lc  N/* 

-?2c  N/» 

30* 

6.26  x  104 

10.99  x  104 

45* 

12.84  x  104 

12.5  x  104 

90# 

25.4  x  104 

Using  the  results  of  Tsble  5  s  plot  of  creek  sngle  (0)  vs.  the  J}  at 
initiation  (Jjc)  was  aide.  This  is  shown  in  Figure  ll.  An  eaplrlcal  fit 
of  these  data  suggested  that  Initiation  occurs  when 


Jj  ■  sin*0  ( ! ) 


Turning  now  to  impact  loading  conditions,  a  dynamic  generalization  of 
Equation  (1)  would  be  given  by 


3,  -  OJ  sln2S  (2) 

1  c 

where  now  contains  inertia  effects  while  D  is  a  rate-dspendent  material 
constant.  If  the  loading  rate  effect  is  small,  then  Dsl,  The  assumption 
will  be  made  here  on  the  basis  that  It  is  consistent  with  the  use  of  the 
static  stress-s train  curve  in  finite  element  calculations.  Therefore, 
strain-rate  effects  are  omitted  from  both  sides  of  Equation  (2).  Now,  using 
Equation  (2)  and  the  results  of  dynamic  finite  element  computations,  a 
prediction  of  the  time  at  which  crack  growth  will  initiate  can  be  made.  A 
comparison  of  these  predictions  with  the  experimentally  measured  time  at 
initiation  is  contained  in  Figure  12.  It  can  be  seen  that  the  prediction  is 
quite  reasonable. 


0 


FIGURE  12.  COMPARISON  OF  PREDICTED  AND  EXPERIMENTALLY  MEASURED 
TIME  AT  INITIATION  IN  IMPACT  LOADING 


DISCUSSION  OF  RESULTS 

The  crack  initiation  criterion  suggested  by  Equation  (1)  is  based  on 
fitting  a  curve  through  only  four  data  points  (three  experimental  results 
together  with  the  origin).  To  exploit  the  possibility  that  some  other 
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function  of  $  could  provide  just  as  good  a  fit,  the  problem  depicted  in 
Figure  13  was  considered.  This  is  an  infinite  plate  containing  a  crack  of 
length  2a  that  is  subjected  to  a  remote  uniaxial  applied  stress  whose 
direction  makes  an  angle  6  with  the  crack  line.  The  linear  elastic  solution 
of  this  problem  gives  [7]: 


Kj  -  o(ira)1^2  sin28 


and 


1  /  2 

KII  "  sin  ®c08  B. 


(3) 


The  relation  between  and  the  stress  intensity  factors  Kj  and  Kji 
will  give 


o(’ra)1^2|sin26 


(4) 


where  y  is  the  shear  modulus,  and  k  *  3-4v  for  plane  stress  and 
k  *  (3-v)/(l+v)  for  plane  strain. 


The  expression  within  the  brackets  is  the  energy  release  rate  G  which,  in  a 
linear  elastic  problem,  la  Just  J.  Therefore,  Equation  (4)  can  be  written  as 


Jj  -  J  sln2B 


(5) 


Now  if  fracture  occurs  when  J  achieves  a  critical  value  that  is  Independent 
of  8»  then  Equation  (5)  will  be  identical  to  Equation  (l). 

Despite  their  similarity,  there  is  a  significant  difference  between  the 
origins  of  Equations  (1)  and  (5).  Specifically,  Equation  (5)  contains  J} 
found  by  taking  the  material  to  be  linear  elastic  and  assuming  a  fracture 
criterion;  l.e.,  J  ■  for  all  B  values.  In  contrast,  Equation  (l)  was 
obtained  by  computing  J j  by  elastic-plastic  finite  element  computations  and 
using  actual  fracture  data.  Two  conclusions  can  tentatively  be  drawn. 

First,  the  functional  relationship  between  ^  and  J  in  the  linear  elastic 
case  may  remain  unaffected  by  plasticity.  Second,  this  criterion  may  be 
geometry-independent.  If  indeed  this  is  the  case,  the  problem  of  predicting 
crack  initiation  reduces  simply  to  knowing  J^c  for  the  material  and  per¬ 
forming  a  linear  elastic  analysis  of  the  problem  to  determine  the  relation¬ 
ship  between  J  and  Jj.  These  implications  present  intriguing  possibilities 
which  saist  be  further  explored. 
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FIGURE  13.  INFINITE  PLATE  WITH 


INCLINED  CRACK 

In  regard  to  the  prediction  of  crack  Initiation  in  impact  loading , 
Figure  12  shows  quite  reasonable  agreement.  However,  the  significance  of 
this  observation  is  not  certain  at  this  time.  The  fact  that  the  predicted 
time  at  Initiation  in  all  three  cases  lies  within  the  measured  time  Interval 
from  the  Instant  when  the  crack  Initiation  was  detected  below  the  surface  of 
the  specimen  to  the  Instant  when  it  appeared  at  the  surface  Is  very  encourag¬ 
ing.  That  Is,  since  the  two-dimensional  finite  element  analyses  used  In  the 
prediction  do  not  model  crack  tunnelling,  the  predicted  time  should  be  com¬ 
parable  to  some  Instant  Inside  the  measured  interval.  While  It  can  there¬ 
fore  be  claimed  that  the  experiments  do  provide  data  that  are  comparable  to 
the  analysis  predictions,  more  experiments  are  needed  to  draw  any  decisive 
conclusions. 

CONCLUSIONS 

The  results  described  In  this  paper  Indicate  that  the  3  parameter  can 
be  used  as  an  effective  criterion  for  crack  initiation  In  elastic-plastic 
mixed  mode  loading  conditions.  Because  a  reasonable  prediction  was  made  for 
Impact  loading  using  quasi-s tat  leal ly-based  values,  It  appears  that  the 
effect  of  loading  rate  may  not  be  too  Important  in  HY-100  steel.  Finally, 
the  results  suggest  the  Intriguing  possibility  that  the  geometric  Influences 
on  mixed  mode  loading  under  elastic-plastic  conditions  may  be  Identical  to 
those  arising  In  linear  elastic  conditions. 
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